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Abstract. We consider the Schrodinger operators on zigzag nanoribbons (quasi-lD tight- 
binding models) in external magnetic fields and an electric potential V. The magnetic field is 
perpendicular to the plane of the ribbon and the electric field is perpendicular to the axis of 
the nanoribbon and the magnetic field. If the magnetic and electric fields are absent, then the 
spectrum of the Schrodinger (Laplace) operator consists of two non-flat bands and one flat 
band (an eigenvalue with infinite multiplicity) between them. If we switch on the magnetic 
field, then the spectrum of the magnetic Schrodinger operator consists of some non-flat bands 
and one flat band between them. Thus the magnetic field changes the continuous spectrum 
but does not the flat band. If we switch on a weak electric potential V — > 0, then there are 
two cases: (1) the flat band splits into the small spectral band. We determine the asymptotics 
of the spectral bands for small fields. (2) the unperturbed flat band remains the flat band. 
■ We describe all potentials when the unperturbed flat band remains the flat band and when 

one splits into the small band of the continuous spectrum. Moreover, we solve inverse spectral 
problems for small potentials. 

1. Introduction 

There are a lot of papers about the electronic structure of carbon materials such as carbon 
nanoribbons, nanotubes and fullerenes during the past two decades because of fundamental 
scientific interest in nanomaterials and because of their versatile electronic properties that are 
expected to be important for future nanoelectronics [DDE], [Haj . [SDDj . Among the carbon 
nanostructures, a simple variation of graphene, ribbons has been studied extensively. The 
graphene nanoribbons with varying widths can be realized either by cutting [Hi] mechanically 
exfoliated graphenes [Nolj . or by patterning epitaxially grown graphenes [ZTSKj. 

In our paper we concentrate on carbon nanoribbons, the recently discovered two-dimensional 
carbon crystal |Nolj . |No2] . The graphene forms a hexagonal graph (the honeycomb lattice) 
embedded in M 2 with all edges of constant length, see Fig. 1-3. In physics one considers only 
two type of ribbons: zigzag and armchair. When a single graphite layer is terminated by 
zigzag edges on both sides, then this layer is a zigzag graphene nanoribbon (ZGNR), see 
Fig.l. When a single graphite layer is terminated by armchair edges on both sides, then this 
layer is a armchair graphene nanoribbon, see Fig.l. Following conventional notation, a ZGNR 
is specified by the number N of zigzag chains along the ribbon forming the width and is 
referred as a N-ZGNR. For example, in Fig.l we have the 4-ZGNR. 

The main interest in physics is to study electronic structure of ribbons, see [DDEJ, [Haj, 
[SDDJ and references therein. Moreover, there are numerous physical papers devoted to 
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graphene nanotubes and nanoribons in an external field, see [KLYHj . [N] , |No3j . |SCLj . [SDDJ 
and references therein. 

There are some mathematical papers about graphene- media, see [KL] . |KL1] . jK], |KuP] 
and |Paj . devoted to continuous models, see also |KS] . |RR] about other models. But in the 
physical literature the most commonly used model is the tight-binding model, corresponding to 
discrete Hamiltonians. The zigzag and armchair nanotubes in external electric and magnetic 
fields (discrete models) were considered in |KKu4j . The zigzag nanoribbons (discrete models) 
in external electric field were considered in |KKu3j . 

We consider now the N-ZGNR, N ^ 1 in external fields, see Fig. 2. In fact we have 
a 2D problem. Firstly, there is a 2D nanoribbon C M 2 . The ZGNRs are assumed to be 
infinite along the x direction (horizontal). In general, the boundary of the ZGNR can be very 
complicated. We consider only the ZGNRs with the simple boundary, see Fig. 2. Secondly, 
an external 2-dim transverse electric field S ex t G M 2 is applied across the ZGNR and along the 
y direction (vertical), see Fig. 2 with the 3-ZGNR in an external electric potential. Moreover, 
the homogeneous magnetic field is perpendicular to the ribbon and along the z direction. 

We will show that the spectrum of the Laplacian (an unperturbed operator) on the ZNRB 
consists of two spectral bands separated by a gap and exactly one eigenvalue with infinite 
multiplicity (a flat band) in this gap. If we switch on an external magnetic field, then we will 
show that the unperturbed spectral bands change, but the flat band (unperturbed eigenvalue) 
will be still a flat band. If additionally we switch on an small external electric field, then we 
will show that the spectral bands slightly change and the flat band will splits into either a new 
spectral band or will be still a flat band. Our main goal is to describe this transformation. 
For applications of our model see ref . in |Ha] , |SDD] . The second part of our paper is devoted 
to the inverse problems. 

Denote a N-ZGNR by T = T^. We consider the ribbon r C K 2 as graph. This graph is just 
a discrete set of vertices V = {v K : x = (n, k) G Z x N p }, where p = 2N + 1, N/v = {1, .., N} 
and a set of undirected edges E such that [v, v'] G E if the vertices v and v' are connected by 
a edge. Each vertex inside the ribbon is connected with some other 3 vertices, the vertex on 
the boundary is connected with 1 or 2 other vertices. If the vertices v and v' are connected 
by an edge, we denote this by v ~ v'. Introduce the discrete Hilbert space i 2 (T) consisting of 
functions / = (/JxeZxNp on the set of vertices V equipped with the norm H/H^n = Yl \ f*\ 2 - 
We define the magnetic Hamiltonian A& on the nanoribbon T in an external magnetic fields 
# = .8(0,0,1) G R 3 by 

(A b f) x = f*>, x=(n,k)eZxN p , p = 2N + 1, & = 
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Fig 2. A horizontal zigzag nanoribbon at N = 3, k G {1, 2, .., 7} in the constant vertical electric field. 



where / = (/ x )zezxN p £ ^ 2 (r) satisfies the Dirichlet boundary conditions 

/n,o = fn, P = 0, n6l (1.1) 

The factor e 111 "-^' is associated with the magnetic field SB and if SB — 0, then e* *'*' = 1. 

We define the discrete Hamiltonian = A& + V on the nanoribbon T (a tight-binding 
model of single- wall nanoribbons, see [jSFJD], [N]) in an external electric potential V and the 
uniform magnetic field SB = B(0, 0, 1) G M 3 . The magnetic field is perpendicular to the plane 
of the ribbon and the electric field is perpendicular to the axis of the nanoribbon and the 
magnetic field. 

Our model nanoribbon T is a graph, which is a set of vertices v nj * and edges E njkt j given by 

v n , 2fc+ i = (VS(2n + k),3k), k G ~N° N 
v n ,2k = {V3{2n + k),3k-2), fceN/ 

where N fc = {1, .., k} C N and N° k = N k U {0}, see Fig. 2 and 3 for the case N = 3. 

Our Schrodinger operator has the form — A& + V, b — G T = R/(27rZ), where the 
magnetic operator A& (see more in Section 4) and the electric potential V are given by 

(AJ) n , 2fc+1 = e- ib ^f^ 2k + e^ n+2 ^f n ^ 2k+2 + e^- fc Vn, 2fe+2 , k G 

/n,0 = fn, P +l = 0, p = 2JV + l, 

(A b f) n , 2k = e-W^+Vf^-t + e- ib( - n+2k -^f n+1>2k ^ + e ib ( 2 " +fc )/ n , 2fc+1 , fc G Njy, (1.3) 

{Vf) H = v k f„, x=(n,k)EZx N p , (1.4) 
where rc G Z and / = (/^) x eZxN p e £ 2 (r) and 

^ — (^fc)? ^ I n fact we consider the 

Schrodinger operator Hf, on the ribbon F (on the set Z x N p ) with the Dirichlet boundary 
conditions f nt o = f n , P +i = 0. Note that (11. 3p gives Hb +2n = Hb for all b G 1. Our electric 



En,k,l — [ v n, 2 fc; v n, 2 fc+l] 

-En,fc,2 = \y n ,2ki V n ,2fc-l] ? n G Z, (1.2) 

E n ,k,3 — [ v n, 2 fcj v n+l, 2 fe-l] 
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potential V is given by v G W, since the electric field is perpendicular to the axis of the 
nanoribbon. We formulate our preliminary result. 

Theorem 1.1. i) The operator H b = A b + V is unitarily equivalent to the operator Jj® 27r ) Jt^, 
where J t = J t {b,v) is a Jacobi operator, acting on C p and given by 



[JtV)n = a n -\Vn-\ + anUn+1 + VnVn, = 

p = 2N + l. 



even n 



2|cos(f 



3n-2 



6)| odd n 



;i.5) 



where y = (y n ){ G C p , y = = y p+1 
ii) The spectrum of H b is given by 

a{H b ) = |J <r k , a k = Ajfe([0,27r]) 

k& N 




[A+A^L k<0 



k G Zjy = {—N, N}, (1.6) 



where A_7v(t) ^ \-N+i(t) ^ ••• ^ Ajy(t) are eigenvalues of J t ,t G [0, 27r). Moreover, A„(-) «s 
rea/ analytic for any t, where a2k-i ^ 0, k G Nat. 

Remark. 1) Below we will sometimes write Xk(t,b,v), aj-(b,v), J t (b,v),.. instead of Xk(t), 
ak, Jt,--, when several potentials v or magnetic fields are being dealt with. The operator 
J t = Jt{b, v) is the symmetric p x p matrix given by 
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2) We take the specific boundary of the ZGNR, see Fig. 2. If we change the boundary of the 
ZGNR, then the corresponding operator J t will be more complicated. In the last case even 
the exact calculation of the spectrum for the unperturbed operator A can be the problem. 

3) If Afc(t) = const for all t G [0, 2tt], then ak is a flat band. Otherwise, ak is a non-flat spectral 
band. In the periodic spectral theory the basic problem is to describe all flat bands. Examples 
of the dispersive curves for Afc(-) are given in Fig. 5. 

We recall the result about the spectrum of a(A ) of the Laplacian A with v = and 6 = 
from [KKu3j (for more details see fl2.3l) - fl2.5p ) 



ct(A ) = a ac (A ) U a pp (A ), a pp (A ) = {0}, a ac (A ) = [-A u , A u ] \ [- Sl , Sl 



where A = (5 + 4ci 



an 



an 



cos ■ 



sin 



a G 



N+r N+r 

Recall that = {— N, N}. Now we describe the spectrum of the magnetic operator A&. 
Theorem 1.2. i) Let b G T. Then the spectrum of At is given by 

a(A b ) =a ac (A b )Ua pp (A b ), a ac (A b ) = (J A fc ([0,27r]), a pp (A b ) = {0}, (1.9) 

kei, N \{o} 

where ao(A b ) = {0} is a flat band and A^(-) is defined in Theorem \l.l\ 
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ii) Let b G T and b — > 0. Then a ac (A b ) = , ji 2 ] U [/i^, /i^] and /i^ nas asymptotics 

Hf = ±(5 + 4d)* + 0(6 2 ), ^ = ± Sl T 2(^? 6 + °^)- ( L1 °) 
We describe the spectrum of H b for b and t> = (i> n )i £ ^ p for which a pp (H b ) ^ 0. 

Theorem 1.3. Let (b,v) G T x W. Then a pp (H b ) ^ iff v 2n +i = v l for any n G Njv. 
Moreover, if a pp (H b ) ^ ; then a pp (H b ) = cr (H b ) = {vi} and each \o- n (H b )\ > and a n (H b ) C 
(v i, +oo), o-- n (H b ) C (— oo, f i) for all n G N/v. 

Remark. (1) Theorem 11.31 generalizes results from |KKu3j . devoted to the operator Ao + V. 

(2) If a pp (H b ) = {vi}, then a flat band o~o(H b ) = {vi} lies in the gap of the continuous 
spectrum, see Fig 3 and 4. 
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Fig 3. Spectral bands for the case v, b = and N = 15. 

^ unperturbed potential v = 

-3 CT_ 2 cr_i CT o cr 2 CT 3 (<jq is a flat band) 



— small perturbation v 

Fig 4. Spectral bands for the case N = 3. 

We will use notations Xk(t,b,v) = Xk(t) for our spectral curves Xk(t). The perturbation 
theory gives the standard asymptotics 

X k {t,b,v) = X k {t,b,0) + O{\\v\\) asv^O. (1.11) 

In Theorem 11.41 we determine asymptotics of the eigenvalue Xo{t,b,v) as v — > at fixed 6, 
recall that Ao(t, o, 0) = 0. Below we will use this result to solve inverse spectral problem for 
small potentials. 

Theorem 1.4. Let b G T and let v -)• 0. Then 

N 

Ao(0 = $> 2fc+ i%(£) + 0(|M| 2 ), (1.12) 

n=0 
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1) unperturbed potential 6, v = 2) small perturbation b, v 

Xo(t) = is a flat band Aq(£) ^ const covers band 



Fig 5. Spectral curves X k (t) and bands a k for the case N = 3. 

Strong electric fields. We consider now the nanoribbon in strong electric fields tV as the 
coupling constant r — > oo and V is fixed. Our operator has the form Hb(r) = + tV . We 
determine asymptotics of the spectral bands a k (jv) = [A£ (r), A^(r)], /c e Zjv of the operator 
Hb(r) as the coupling constant r — > oo. 

<ti q j 4 05 q_6 

TV\ TV2 TV3 TV4 TW5 Tt>6 A 

Fig. 6. Spectral bands for the case of strong electric field 



Theorem 1.5. Let H h {r) = A b + rV , where r — > 00 and let the vector v = (ffc)i G R p ,p = 
2iV + 1 satisfy v± < .. < v p . Then the spectral bands o-j(rv) = [A^(r), A~(r)],j G Z N have 
asymptotics: 

A- (Tj=TV fe , 4 fc = 1 , s = {-I) , k = j + N + 1, 



r v k _i - v k v k+ i - v k 

p+1 



t; = t; p+1 = 0, rf = r±=0, r± +1 = 1, r 2n = 0, r+ = 4, nG%, (1.14) 



kjW|= 4 + Q(r " 1) ^A/, anrf KHhO(r- 2 ). (1.15) 
Moreover, the operator i?&(r) /ias not a /?a£ 6an<is /or r Zarye enough. 
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Remark. 1) If r — > oo, then roughly speaking the spectrum of operator Hj,(t) consists of 
p bands o~k{Tv),k G N p with lengths \ak(jv)\ > 0, separated by 2N large gaps. In this case 
all 2N gaps are realized. In same sense it is clear, since we have the operator — = V + eA& 
with small coupling constant e = -. The operator V acting on £ 2 (T) has only flat bands 
(eigenvalues) {vk},k G N p . Under the small perturbation eAb these flat bands {vk},k ^ N 
become the small spectral bands ak(rv) of the continuous spectrum with lengths > 0. It is 
important that only Theorem 11.31 shows that apf(rv) is not a flat band, since asymptotics 
(jl,15p are not sharp for j = N. 

Inverse spectral problem for odd potentials. Below we solve inverse problems for 
sufficiently small "odd" potentials. Define the space of odd potentials 

V odd = {v = {v k )l e R p : v 2k = 0, k e N N } 

equipped with the standard norm || • || in W. Let (b, v) G T x V dd and let ^ to < ti < ... < 
^ 7r. We define the mapping A : V m ~ > 7& N+1 by 

v i-> A(v) = (X (tj, 6, v ))q , where \ (tj,b,v) E a = A o ([0, 2ir], b, v), 

and Ao is defined in Theorem 11.11 Let B r — {v G V dd '■ \\v\\ < r}, r > be the ball in V dd- 

Theorem 1.6. i) Let b G T be sufficiently small. Then the mapping A : B r — > A(B r ) is a real 
analytic bisection for some r > 0, where A(B r ) C R^ 1 is an open domain and G A(B r ). 

ii) In the case of i) for any e G (0, r) there exists a potential v G B £ and a potential w £ B £ 
such that A (it?) = A(v), i.e. there is no a global injection. 

Hi) Let b G T and let tk G [0, 27r), k G N p be a sequence of a distinct numbers. Then the 
vector (Ao(tfc, b, v))i uniquely determines the potential v G B r for sufficiently small r > 0. 



A 




Fig 6. The case N = 3, small perturbation b, v. Spectral curve Xo(t) covers the band do; 
2(N + l)-periodic eigenvalues Afc = Xo(j^j) belong to cfq. 

Remark. Let t k = k G N° N . Then X (t k ) are so-called 2(A^+l)-periodic eigenvalues for 
the Jacobi matrix J t . Thus, if b, r are small enough, then by Theorem II .61 i). 2(A^ + l)-periodic 
eigenvalues Ao(t&) G ao, k G N% uniquely determine the odd potential v (and operator V). 
Note that we take only some part of 2(A^ + l)-periodic eigenvalues to determine the potential 
V, since Xj{tk), j ^ are 2(A^ + l)-periodic eigenvalues too. 
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Inverse problems for monotonic electric potentials without the magnetic field. 

We consider the case 6 = 0. Define the mapping in terms of antiperiodic eigenvalues 

*:T,^*(T,) = (A fc (7r,0, 

Here Afc(vr,0,f) is a antiperiodic eigenvalues for Jacobi operator J t , since t — ir. Define the 
bounded set of monotonic potentials by 

V a = {v G M. p : ^ Vi < v 2 < v 3 < .. < v 2 n+i ^ «}, 

V a = {v G W p : a ^ vi > v 2 > v 3 > .. > v 2 n+i ^ 0}. 

Theorem 1.7. i) The mapping : V a — > M> N+1 is an injection iff a G [0, 1]. 

ii) The mapping ^ : V a — t ^ N+l is an injection iff a G [0, 1]. 

Hi) In these both cases there exist the algorithm of recovering of v in terms of $f(v) = 

iv) Spectral curves Xk, k G Z^v satisfy Xk(7i+t, 0, v) = Xk(7r—t, 0, v) for any (t, v) G [0, n] xf p . 
If v G Vi or v G Vi then all A / t(7r,0,f) ; k G Zjv are distinct numbers. 

Remark. 1) This Theorem shows that the set of all antiperiodic eigenvalues uniquely 
determine the increasing (or decreasing) bounded potential. Moreover, in the proof of Theorem 
11.71 we show how to determine the monotonic potential v by antiperiodic eigenvalues ty(v) = 
(X k (ir,0,v))^ N . 

2) Theorem ll.7l iv) shows us that reach the maximum or minimum at the point t = n. 
Then for sufficiently small monotonic potentials v the spectrum change multiplicity at any of 
anti-periodic eigenvalues Afe(7r, 0,v). 

3) Note that in the paper [KKu2] we proved that some finite set of fc-periodic eigenvalues 
uniquely determines the spectrum (including multiplicity) of a periodic matrix-valued Jacobi 
operator. In the present paper Theorem II .6l i gives the stronger result and shows that some of 
iV + 1-periodic eigenvalues uniquely determine the operator Hb for potentials v small enough. 

In the proof of our theorems we determine various asymptotics for periodic Jacobi operators 
with specific coefficients given by (11.51) . Note that there exist a lot of papers devoted to 
asymptotics and estimates for periodic Jacobi operators, see e.g. [KKulj . [KKrj . |Laj . |vMouj . 

We present the plan of our paper. In Section 4 we prove Theorem 11.11 with a technical proof, 
where we use arguments from |KKu3j . In Sect. 2 we describe the properties and prove 
Theorem ll.2[ In Sect. 3 we prove Theorems 11.31 - fTTTl 



2. Properties of the magnetic operator A b 

Theorem 1.1 is proved in Section 4. In this section we describe the spectral properties of 
A;,. Recall that Theorem 1.1 gives that the operator = A5 + V is unitarely equivalent to 
the operator Jj^s ^||, where J t = J t (b,v) is a p x p Jacobi operator. Introduce a Jacobi 

operator J t = J t (b,v) acting in C p and given by 

— 1 even n 

(Jty)n = a n -xVn-\ + a n y n+1 + v n y n , a n = < 3„- 2t \ > (2- 1 ) 

I 2 cos(| — m 2 - ± b) odd n 

and recall that J t = J t (b,v) is given by 

(Jty)n = a n _!y n _i + a n y n+1 + v n y n , a n = \a n \, (2.2) 
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where y = (y n )i £ C p and yo = = y p +i,p = 2N + 1. Note that the matrices J t = Jt(b,v) 
and J t = Jt(b,v) are unitarely equivalent. Then the eigenvalues \k(t) = Xk(t,b,v), k G 7L^ 
of Jt (see Theorem ll.ll ii) are eigenvalues of Jt ■ Sometimes below we will use the matrix J t 
instead of J t to study the functions Xk(t) and the bands a k , defined in Theorem ll.ll ii). 

The Laplacian A . We describe the spectrum cr(A ) of the unperturbed operator A 
at V = and 6 = 0. Let X°_ N (t) ^ X°_ N+1 (t) ^ ... ^ X%(i) be eigenvalues of the Jacobi 
matrix J 4 = J t (0, 0) corresponding to A (see Theorem II. H i). We recall the results about the 
spectrum of a(A ) from |KKu3] : 

i k t 
X° k {t) = (a 2 - 2c k a + 1)3 — , 1 < |A;| ^ N, X° (t) = 0, tG[0,2vr], a = 2|cos-|, (2.3) 

\k\ 2 



N 

' ac w " pp 



*° = at U at = |J at a% = a° Q = {0}, a° ac = [-X° N (0), X%(0)\ \ [s h 8l ], (2.4) 



-N 



*-fe = -0*> ^fc= A fc ,A fc = < , allk = l,..,N, (2.5) 

[[Afc(7r),A^(0)], z/c fe ^0 

J _ rx o,- x o I+1 n /[A fc (t fc ),A fc (0)] = [ Sfc ,(5-4c fc )^], if c k >0 

k ~ i k ' k 1 ~ - k ~ [AM, Ag(0)] = [1, (5 - 4c,)§], i/ c fc < ' (2 ' 6) 



for all k = 1,..,N, where X° k (t° k ) = s k , 2cos| = c fc , and A£(?r) = 1, A^(0) = (5 + 4c x )i 

We use notation J t (b,v), J t (b,v) and Xk(t,b,v) for matrices J t , J t (defined in (12.11) . (11.51) ) 
and its eigenvalues. 
For k 7^ introduce 

*l{t) = (^,n(0)L 1; = ( ^^ nfc ' ^.*»(* + *) = VN + l- k ,2n(t), * G (0, 7r), 

vU-i(*) = ( " 1)W+1 ^^u|^ (n " 1)fc) - ^-x(* + *) = ^i-^xW, * e (0, VT). 

(2.7) 

Below we use notation (■, •) for the standard scalar product in C n . We consider the case V = 
and b — > 0. 

Lemma 2.1. z'j The following identities hold true 

j?m)=^(mi(t), ih&2(f)ii = i, k^o, te[oM, (2.8) 

i.e. & k (t), k are orthonormal eigenvectors for J t °. 
zz ^ For any t G [0, 27r] \ {vr} iae following identity holds true 

sin - / ft \ t\ 

d b (X k )(t, 0, 0) = (jV + 1) 2 A o (t) ( + l)(3iV + 1) (2 cos - + c k sign(t - vr) J - 6 cos - . (2.9) 

Proof. Direct calculations gives i). 
ii) Using the perturbation theory, we obtain 

t - 

(X k y b (t,0,0) = <(Ji)i(0,0)*2(t),*2(t)> = 2sin-]T(6n-5)< 2n _ 1 (t)< 2n (t) 

n=l 
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2 sin § „ t 

-q— 2^(6™ - 5)(2s nfc cos - - s nk s {n _ l)k ) 



2 sin 



A' 



(N + 1) A°(i) Y^ 12ns ^ cos 2 ~ 10s ™ fc cos 2 ~ 6ns « feS («-i)fc + 5s «fc s (n-i)fc) 



2 sin | 



(N + l)\° k (t) 
where we use identities 

N N 

C 2n fc+a = S 2n fc 



(iV + 1)(3JV - 2) cos ~ - (N + 1)(3JV + l)y 



1 < |Jfe| < iV, ael, 



n=0 



n=0 
N 



N 



1 - c 2nfc iV + 1 



n=l 



n=l 



iV + 1 



"Cfe, 



n=l 
AT 

22 nC2kn 

n=l 



n=l 
iV + 1 

A? 



iV 



2J nS 2fcn 



n=l 
iV 



iV+ l c fe 
2 s k 



1 < |Jfe| < N, 



n- 



l-c 2nk (N + l) 



n=l 



n=l 



iV 



N 



N 



n=l n=l n=l 

(JV + l)(JV + 2) 
= i Cfc. ■ 

Lemma 2.2. Lei b,t E [0, 27r] and f 2n +i = /or a// n G N^-. T/ien A = is a simple 
eigenvalue of Jt- 

Proof. Define a Jacobi matrix A by 

/ ax .. \ 

ai v 2 a 2 



A 



a 2 








V 



a 2k 7^ 0, for all k G N n , 



.. a 2n / 

for some n G N. Then A = is a simple eigenvalue of A and x = (x k )1 n+1 is a corresponding 
eigenvector, where x\ = 1, x 2k+ i = (— l) fc ai " ct2fc ~ 1 , x 2 fc = 0, G N n . ■ 

Proof of Theorem 11.21 i) Lemma 12.21 gives that {0} is a flat band. Then Theorem 11.31 
i) and the Theorem 11.11 give the proof of i). 

ii) Relations (1Q]) - (I23]1 imply 

/4(0) = min AS(f) = A°(*°) = A°(z 2 °) = Sl < A°(i), f G [0, 2tt] \ {z x V 2 }, (2.10) 

iG[0,27r] 
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where z\,z\ G (0,27r), z\ ^ z\ such that 2|cos^-| = ci, j = 1,2. We denote A°(t,6) = 
Ai(t, 6, 0), Perturbation Theory gives us that X°(t,b) is an analytic function in some neigh- 
bourhood of points t = z®, z 2 and 6 = 0, because A°(^i, 0) and Ai(z2, 0) are simple eigenvalues 
of analytic matrix function Jt(b, 0). Using 

ftA?(4,0) = 0, #\?(*J,0)>0, j = 1,2, (2.11) 

and the Implicit Function Theorem we deduce that there exist only two functions Zi(b), z 2 (b) 
(analytic in |6| ^ 5 for some 5 > 0) such that 

zj(0) = z° p d t \° 1 (z j (b),b)=0, d*\° 1 (z 3 (b),b)>0, j = 1,2 (2.12) 

for all \b\ 5. Identities (I2.10I) - (I2.12I) yield that Zj(b) are points of local minimum of function 
A?(t, 6) and one of these is a point of global minimum, that is 

u£(b)= min A?(t,6) = min{A?(z,-(&), &)}, (2.13) 

A?(^(6),6) = A?(zJ)+6(9 6 A?(^(6),6)| fe=0 )+O(6 2 ) = Sl +6(«9 6 A?(z#), 6)| 6=0 ) +0(6 2 ), (2.14) 

^A?(^(6),6)| b=0 = d t \1(zl 0)d bZj (b) + $A?(^,0) = &A?(z?,0), j = 1,2. (2.15) 

Then, substituting d 6 A?0°, 0) = <9 fe Ai(z°, 0, 0) from (F23]) to (I27T51) . substituting fl^T3]) to (I2TTD 
and using (12.131) we obtain asymptotics for fi 2 - Other cases are proved similar. ■ 

3. Proof of the main results 



In order to prove Theorem 11.31 we need two Lemmas. 

Recall that v = (ffc)i, p = 2N + 1 and a 2k — 1, k G and a 2k +i = a,2k+i(t,b) = 
2 cos(| — A; G Njy. These functions are analytic for t G C, below we will consider some 

cases when t G C not only £ G [0, 27r] as we did above. Introduce the matrices 

T h (t,b,v)= (~ l ~ V2k ~~2 \ keN N+1 , v 2N+2 = 0. (3.1) 

Introduce the functions Uk = Uk(t,b,v), k G N° N+1 by 

Uo ) = ( ? Y f U2fe W( U2fe - 2 Y fcefW (3.2) 

Note that is an analytic functions for t G C. Below we need a simple fact. 

Lemma 3.1. Lei £, 6 G [0, 2ir] and let v G M p . Suppose a 2 k~i(t,b) ^ /or a// A; = 1,..,N. 
Suppose A = is eigenvalue for J t (b, v). Then u 2 N+2(t, b,v) = 0. 

Proof. Let Jt(b, v)F = XF = 0, where F = (f n )i 7^ is the corresponding eigenvector. Then 
the definition (12.21) implies 

vifi + a 2 f 2 = 0, a 2 fe_i/ 2fc _i + v 2k f 2k + a 2k f 2k+1 = 0, k G Nat, a 2N f 2N + v 2N+1 f 2N+1 = 0, 

which yields 

h = Mi, f / 2fc ) = J-T*f ) = f f ? ; 2fc Y k g Njv+i, 
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6c = rin=i a 2n-i which yields fi 7^ 0, since Then 

n _ ~ r r _ a 2NU2N + V 2 N+\U 2 N+1 _ U2N+2 

U — a 2 N J2N + u 2jv+i/2jv+i — 77 — — 77 — i 

JlQN JlsN 

which yields U2N+2 = 0. ■ 

Lemma 3.2. Let b E T, v E M p . T7ie following asymptotics hold true 



'JV+lj 



|«2fc+i(**)| =e 2H + o(e 2 <*- 1 > t Y t^+oo, z = V C l, fcGNiv. (3.3) 

Proof. Note that 

a k (it) = e i + 0(l), t^+oo, G N p . (3.4) 
For k = 1 identities (E3D and (j3TT]) -(l3T2 ]) yield 

Wit) J - 1 H ^ ) V 0(e( fe -^) J' t ^ + °°- (3 ' 5) 

Substituting ( 13 .5p for fc = 1 into (13. 2p and using (13. ip . ( 13. 4 p we deduce that (I3.5P is true for 
k = 2. Repeating this procedure we deduce that (13. 5p is true for any k G N p , which yield 
asimptotic (13. 3p . ■ 

Proof of Theorem 11.31 If v\ = then Sufficiency follows from Lemma 12.21 If v 1 7^ then 
the proof is similar, since we may consider new operator Hf, = + (V — v\I), where I is an 
identity operator. 

Necessity. We have a pp 7^ 0. Without lost of generality we may assume G cr pp (in other 
case we shift spectrum by adding some constant c to any of diagonal components v^). Then 
A = is an eigenvalue for matrix J tip, v) for tE.Ii, for some set I± C [0, 2tt] and > 0. There 
exists infinite set I 2 C h, satisfying Ylr^i^n-iit, b) 7^ for any t G I 2 . Then, by Lemma l3Tl] 
we have u 2 ^ +2 {t,b,v) — 0, t E I 2 . Then u 2 ]y +2 {t,b,v) = for any t E C, since u 2 n +2 is an 
analytic function by t. Using (I3.3P we deduce |m2at+2(^)| ~ |^2jv+i|e 2Ar< , t — > +00, which yields 
V2N+1 = 0. Then, using (13. ip . (13.21) and f 2 Ar+i = 0, we obtain u 2N (t, 6, v) = -u 2N+2 (t, b, v) = 0, 
t E C. Similarly to above, we deduce v 2 n-i = and so on. 

By O we know that A_ fc (t, 0, 0) < X (t, 0, 0) = < X k (t, 0, 0) for any (t, k) G [0, 2tt] x N n . 
Let some b G T and let v E W be such that v 1 = V3 = .. = v 2 n+i = 0. We fix t E [0, 2n] and 
let variable r G [0, 1], perturbation theory gives us that Xk(t,rb,rv), k E 1*n are continuous 
functions depending on r G [0,1], since these are eigenvalues of analytic matrix-function 
J t (rb,rv). Note that X (rb, tv) = for any r G [0,1]. Then inequalities X-k(t, rb, rv) < 
\ (t,rb,Tv) = < Xk(t,rb,Tv), k E Njv remains for any r G [0,1], since by Lemma I2T21 
A = is always simple eigenvalue of J t {rb,rv). Then by (11.61) a_ k {b,v) C (— 00, 0) and 
o~k{b,v) C (0, +00) for any k E Njv- If V\ = v 3 = .. = v 2 ^+i 7^ then we just shift spectrum 
by adding —v% to any of diagonal components Vk of the matrix J t and repeat the reasoning. 
■ 

Proof of Theorem 11.41 If V = 0, then Theorem 11.21 gives that o~o = {0} is a flat band 
for H b . Also Aq = is a simple eigenvalue of matrix J t (b,0). Using perturbation theory, for 
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v — >• we obtain 

2JV+1 

Xo(t, b, v) = A (t, b, 0) + ( d ^o){t, b, 0)v n + 0{\\v\\ 2 ), uniformly for t G [0, 2tt], 

n=l 

(^Ao)(t,6,0) = ((d Vn J t )(b,0)^ b (t),^(t)} = (^) 2 (t), 

where = (<p„(t))i is a corresponding eigenvector for an eigenvalue Ao(t, b, 0) = 0, i.e. 

Jt(6, 0)$ 6 (t) = and = 1. Using results from the proof of Lemma 12.21 we obtain 

VzrSf) = 0; P^n+iit) = ( — ^) n fin(t) (recall that (3 n are defined in (11.131) ). since J t (b,v) has the 
form (pi~TT|) . ■ 

We start to prove our results about the inverse problems. It is convenient to prove Theorem 
11.61 after the proof of Theorem 11.41 since we use the same notations. 
Proof of Theorem 11.61 We use notation from Theorem 11.41 

i) Let 6 = 0. For any fix t G M the function Xo(t, 0, ■) is real analytic for small arguments, 
since Ao(t,0, 0) is a single eigenvalue of analytic matrix-function J t (0, 0). Then for small 
v G Vodd there exist Frechet derivative d v A(v) : V dd — > M. N+1 is a real analytic function. Using 
([I3ZD , ([USD we obtain 

D = d v A(0) = (Vm(tn))t N 0, m =0 G R( JV "+ 1 ) X ( iV+1 ) ) 

where ^ m are defined in (I1.13p . We need to show that D is an isomorphism, or we need 
to show that D is an injection, since dimVodd = dimR N+1 . Suppose D v° = for some 
v o = (y^P e V. Then, using I02J) . ffl~13|) and we deduce 

- / A 2 

/(t n ) = 0, n G N^, w/iere /(*) = a m t; 2m+1 , a = I 2 cos - J . (3.6) 

m=0 ^ ' 

Then / = and t>2 m +i = 0, m G Njjy-, since this is a polynomial of degree N has A/"+ 1 zeroes. 
We have D is an isomorphism and A(-) is a bijection for small arguments. If b ^ is a 
sufficiently small value then D is an isomorphism too, since D analytically depends on b, 
then A is a bijection for small arguments. 

ii) Define the potential v° = (f°)i G V dd, where v® = — 1, t>2 n+1 = n, v® n = 0, n G Nat. Let 
s > be sufficiently small. Then using (II. 5p we obtain 

J t (b, 8~ V) = s-U(s), A(s) = diag(t; ) + sJ t (b, 0). (3.7) 

Let X-n(s) ^ .. ^ Aat(s) be eigenvalues of the matrix A(s). We have A(0) = diag(t>°) and 
then A_at(0) = -1, A_Ar +n (0) = 0, if n G Nat and A n (0) = n, n G N N . 
Consider the asymptotic for \_ N+n (s), n G for s — » 0. We have e n = (<52nj)j=i (where 
is a Kronecker symbol) are eigenvectors of A(0), corresponding to eigenvalues A_Ar +n (0). 
Using famous results of the perturbation theory, we obtain that the numbers A'__ 7V+n (0), n G 
Nat are eigenvalues of matrix A = ((A' (0)ej ,ek))f tk=v Easy calculations gives A = 0, then 
X'_ N+n (0) = 0,ne N N , and then 

X- N+n (s) = A_jv+n(0) + s\'_ N+n (0) + 0(s 2 ) = 0{s 2 ) ass^O, n G N N . (3.8) 

Using the fact, that X n (t, b, °), n G are eigenvalues of J t (b, sv°) and using (j3.7j) . (13.81) 
we obtain 

A (t, 6, s~V) = s _1 A (s) = O(s), s^0. (3.9) 
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Then A(s _1 f°) = O(s), s — > 0. Using i) we obtain, that for any sufficiently small s > 
there exists unique sufficiently small potential w s G V dd, which satisfy A(s^ 1 v°) = A(w s ), and 
w s s~ 1 v°, since s _1 t>° has large norm. 

hi) Define the mapping A(v) = (Xo(t k , b, f))o- We use similar arguments as in i). For any fix 
(t,b) G M 2 the function Aq(.) = \o(t,b, •) is real analytic for small arguments, since for small 
v the number Xo(t, b, v) is a single eigenvalue of analytic matrix-function J t (b,v). Then there 
exists Frechet derivative d v A(v) : V a dd M. 2N+1 , which is real analytic function for sufficiently 
small v G V dd- Using (I1.12p . HI . 13|) we obtain 

A, ^ a,A(0) = (v. m+ i(tn))llh N m =o e E (2Ar+1)x(JV+1) . 

We need to show that the matrix (operator) D is an injection. Suppose D v° = for some 
w o = (^) 27V+1 G V QfM . Then, using ffTT2|) . (Q3I) . we deduce 

TV 

/(*„) = 0, n G N^, w/iere /(*) = ^ (3.10) 

m=0 

Using (I1.13P and (11.51) we obtain 

m 

= $^( c J> m cos & + s i' m sin ^)' Cm > m 7^ 0, m G (3.11) 
i=o 

for some constants Cj >m , Sj^ m . Then using (13.101) . (13. lip we deduce that f(t) = 0, since / is a 
trigonometrical polynomial of degree N and it has 2N + 1 zeroes in the interval [0, 2tt]. Also 
we deduce that t>2m+i = 0, m G N%, since by (13.111) we have that /3^(t), m G N° N are linearly 
independent functions. Then v = and Do = d v A(0) is an injection. Also A(-) is an injection 
for small arguments, since d v A(-) is a real analytic function for small arguments. ■ 

We consider the case of the strong electric fields. 
Proof of Theorem 11.51 In the case A& + tV the corresponding Jacobi operator depends on 
r and is given by 

(J t {b, Tv)y) n = a n -iy n -! + a n y n+ i + rv n y n , y = (y n ) n& G C p . (3.12) 
Using (12.11) we rewrite J t (6, rv), t — > oo in the form 

J t (b, tv) = t(V + r^JtQ), 0)) = r(V + eJ t {b, 0)), V = diag(u.-)?, as e = - 0. 

T 

Then the perturbation theory (see Sect. XII, 1, [RS] ) for V + sJ t (b, 0) gives 
A fc _AT_i(t, 6, rv) = r(v k - a k e 2 + 0(e 3 )), 

u 2 

a k = V — w fcli = (e; j J t (6,0)eg), (3.13) 
jeN p \{k} J 
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Uk,j = Uk,j{t,b) 



2 2 ")■> j N ' f3 14) 

2cos(|-^6), j = k- 1, k e 2N v ' ; 



for any & G N p , here Ve^ = Vje® and the vector e° = (5j, n )£ = i G C p . The matrix J t (b,0) 
{ukj} is given by ( 12. ip . where 

fO, 

2cos(|-^6), j = k + l, JG2N, 

ot/ier cases 

Then (l37L3"p . (ETiip yield (Oljl . ( Q5j) . ■ 

Proof of Theorem 11.71 i) Sufficiency Sufficiently to show that \I/ is an injection for a = 1 , 
since Vp C V a for (3 ^ a. Using (II. 5p we get a 2n -i(rr,0) = 0, a 2 „ = 1, n G Njv, so direct 
calculations give 

N N 

[j {\ k (7r,0,v)} = a(J 7r (0,v)) = {v 1 }u\J{\: (A - v 2k )(X - v 2k+1 ) - 1 = 0}. (3.15) 

fc=-iV k=l 

Below we use the notation A^ = Afc(7r, 0, v). Firstly, let Pj(X) = (A — Xj)(X — y^) — 1, j = 1, 2 
be two quadratic polynomials for some ^ x± < yi < x 2 < y 2 ^ 1. Let //j < i/j be roots of 
Pj, then it is not difficult to show that 

Mi < 1^2 < < 1 < v\ < u 2 . (3.16) 
Using this fact, monotonicity of v and (13.151) we deduce that 

X% < .. < \ {0 \ < < Ai 0) < 1< \® < .. < A<J>, (3.17) 

where A^ = V\ and A_^_ 1+fe , are roots of the polynomial (A — t>2fc)(A — v 2k +i) — 1 for 
k G Nat. Then ty(v) uniquely determine monotonic potential v and then \l/ : V a — > M? is an 
injection. 

Necessity. Suppose that a > 1. Define two potentials v,w G V Q by 

1 3 

< V! = 2e < - < v 2 < .. < v p - 2 < - < 



Vp-l 



1 + 5e - V2e + e 2 < v p = 1 + 5e + V2e + e 2 < min|a,^|, (3.18) 



1 3 

< w 1 = Ae < - < w 2 = v 2 < .. < w p _ 2 = v p ^ 2 < - < 

Wp _ x = 1 + 4e - 2Ve + e 2 < w p = 1 + 4e + 2Ve + e 2 < minja, ~ j, (3.19) 

for some e > small enough. Direct calculations give us that V(v) = #(iu) = (\ { k 0) ) N N} where 
A_^_ 1+jfc , A^ are roots of the polynomial (A — t>2fc)(A — ^2fc+i) — 1 for A; G ~Nn-i and \®\ = 2e, 
= 4e, AjJ } = 2 + 6e. Then * : V a -> M 2JV+1 is not injection, since = %(w). The 

Proof of ii) is similar to the Proof of i). 

iii) Firstly, Aq '* = v i determine the first component of potential. By the remark after (I3.17P 
the components v 2k , v 2k +i are defined uniquely as a roots of polynomial (A — A^°|f_ 1+fe )(A — 

Af) + 1. 
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iv) By ( II .5p we have that J n+t (0,v) = J n ^ t (0,v), which yields symmetry \ k (ir — t,0,v) = 
\ k (ir + t, 0, v ) and then t = n is a point of local extremum for functions X k . If v G V\ (or V±) 
then by (13. 17j) all components are distinct numbers. ■ 

4. Appendix, proof of Theorem 1.1 

For the magnetic field 38 = 5(0,0, 1) G M 3 the corresponding magnetic vector potential is 
given by 

1 B 

s*{x) = -[38, x] = — (-x 2 ,xi,0), x = (x 1 ,x 2 ,x 3 ) G K 3 . 

Define the coefficients a v ,p{t) = {s^{x v + te U)P ), e u>p ), where u, p G Z x N p and t G [0, 1]. In 
Lemma [4.11 we will show that a u>p (t) does not depend on t G [0, 1]. Recall that the magnetic 
operator A b is given by 

(A 6 /) w = e m "^f n , 2k + e ta ^-^^f n . h2k+2 + e M -^+ 2 >/ n , 2fc+2 , k G N%, 

f f n h- B ^ 

Jnfi — Jn,2N+2 — U, — — - — , 

{A b f) a = e^-W-Vfek-i + e ia ^> 2h -Vf n+1 , 2k -i + e ia ^"^ f n , 2k+1 , k G N N , (4.1) 

where u = (n, 2k + 1), a = (n, 2fc) G Z x N p and / = (f n ,k)(n,k)eZxN p G ^{^)- 

Lemma 4.1. Let a function a u p (t) = {st (x„ + te VtP ), e u>p ), where u, p G Z x M p and t G [0, 1] . 
Denote u = (n, 2k + 1), a = (n, 2A;), i/ien 

a w ,(n,2fc+2) = K n _ V(n-i,2fc+2) = + 2fc), a W) (n )2 fc) = -b(2n + A;) , (4.2) 

a<r,(»i,2fc-i) = ~b{n- k + 1), a CT) ( n+1)2fc _a) = -6(n + 2A; - 1), a CT ,(n,2fe+i) = &(2n + /c), (4.3) 

w /iere b = ^ and all (t, n, k) G [0, 1] x Z x N p . 

Proof. Identity ^(r) = f [e , r], e = (0, 0, 1), r G M 3 yields for any t G [0, 1] 

5 B 



where a = (n, 2k + 2), (n - 1, 2A; + 2), (n, 2Jfe). 
Recall x w = (V3(2n + fc), 3ife, 0). If a = (n, 2A; + 2), then = ±, 0) and gj) yields 

1 x B ( ~2 \ °\ y/ZB . 

V(n,2fc+2) = -w ([e 0) x w J, e W)(T ) = — det 1= — — (n - k). 

\VS(2n + k) 3k 0/ 

If cr = ( n - 1, 2k + 2), then = (-^, §, 0) and (JHj) yields 

ao;,(n-i,2fc+2) = — ([e ,x w J,e Wi(T ) = — det 1 =— — (n + 2fc). 

VV3(2n + ife) 3A; 0/ 



If a = (n, 2A;), then = (0, -1, 0) and (03]) yields 

i 
i 

V3{2n + k) 3k 



B B I _1 \ fiv/3 
a^,(n-i,2fc+2) = — ([e , xJ,e Wj(T ) = — det [ 1|= _(2n + fc). 
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The proof of other cases is similar. ■ 

Substituting identities from Lemma 14.11 into the magnetic operator A& is given by (14. ip we 
obtain (II. 3p . 

Proof of Theorem 11.11 i) Recall that 



(A 6 /) n ,2fe+i = e 

{A b f) n ,2k = e 



-ib(2n+k) 



fn,2k + 6 



ib(n+2k) 



fn,0 — fn,2N+2 ~ 0, 

ib{n-k+l) f , -ib(n+2k-l) r , 

Jn,2k-1 + e Jn+l,2fc-l + e 



2k- 1 



/n-l,2fc+2 + e l6< - n k ' f n ,2k+2 + ^2A;+l/n,2fe+l? & £ N^, 
i6(2n+fc) 



/n, 2fc+l + V2kfn,2k, k G N N , 

(4.5) 



Define the unitary operators [/, 5 acting in £ 2 (Z) by 



r 



Jb 



Sh — (h n+ i) ne z, h — (h n ) 



(4.6) 



For each (f n ,k)(n,k)eZxN p € ^ 2 (H we introduce the function = (f n ,k)n& € ^ 2 (Z), fc G 
= 2iV + 1 and ^ = (^ fc ) feeNp G (£ 2 (Z)f. Using (gSD and (V/) n>fc = v k f n , k for any 
(n, k) G Z x N p , we obtain that the operator H : £ 2 (r) — >■ £ 2 (r) is unitarily equivalent to the 
operator K : £ 2 (Zf -> ^ 2 (Z)*>, given by 

(#VWi = r fc (fT) V 2fc + + 5V 2fc+1 )^ 2fc+2 + v 2k+1 ip 2k+1 , = VVh = 0, k G N^, 

(AT^) 2fc = t/V- 1 + Sf 2 *- 1 )^*-! + U 2 r k ij 2k+U k G N*. 
We rewrite K in the matrix form by 



( v, 


(1 + tS*)U 








\ 




U*(I + tS) 


V2 


U 2 T 








^2 





{U*) 2 T 


V3 


(f+T 3 S*)U 





^3 










U 2 T N 




V 






(f/*) 2 ^ 


v P J 





(4.7) 



where ^ fc G £ 2 (Z). Note that iT = K, since S* — S . We rewrite .fT in the matrix form by 

(vx A 1 .. \ 



AT 



A* w 2 A 2 
A* 2 v 3 A 3 



A 







A* 






p-i 



(4.8) 



/ 



where A 2k = r k U 2 , fceN w and A 2fe+ i = (r k +T 2k+1 S*)U, k G N^-. Define the unitary operator 
% = diag(u fc )i, where u 2k = U 3k ~ 2 , k e N N and u 2k+1 = U 3k , k G N° N . Using US* = rSU* 
we obtain 

(v x r x \ 

i; 2 r 2 .. 


r p -i 



r* w 3 r 3 
. .Or* 



(4.9) 



where r 2k = r k , k G N^v and r 2fe+ i = T- k (r 6k+1 S* + 1), k G N^. 



/ 
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Introduce the unitary operator $ : £ 2 (Z)p -»■ f® J^g, = C p , by $(V>*)i = (<Mfc)i> 
where : £ 2 (Z) — > L 2 (0, 2ir) is an unitary operator given by 

int 

<ph = J2 ; h = {K)n& e £ 2 {Z), t g [0, 2tt]. 

Then we deduce that 



2tt 



[0,2tt) 



dt 

; 2? 



where the operator J t : — >■ has the matrix given by 









1 + Te~ u 





.. 


\ 


(yi\ 






1 + re lt 


t'2 


T 


.. 




V2 


= 







T 


V3 


.. 




V3 




V 










r N v p 


J 


\y P J 



y = {y k )l G C. 



(4.10) 



The matrix J t is unitarely equivalent to the matrix J t , given by 





( 




ai(t) 





.. \ 








a x {t) 




a 2 (t) 







Jt = 







a 2 (t) 
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where a 2fc (t) = 1, keN N and a 2fe _i(t) = 2| cos(| - 
Thus we deduce that the operator H is unitarily equivalent to the operator J^ 27T - } Jt^- 
ii) Let Xf.(t),k G Z N be eigenvalues of J t satisfy A_7v(t) ^ A_Ar +1 (t) ^ ... ^ Ajv(£)- From 

the spectral theory of Jacobi operators [vMj we have that if all af.{t) ^ for some t, then 

A_Ar(i) < A_Ar + i(t) < ... < Ajv(£) and perturbation theory gives us that any of these functions 

is analytic in some neighborhood of t. m 
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